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Context. We present 31.2 days of nearly continuous MOST photometry of the roAp star 10 Aql. 

Aims. The goal was to provide an unambiguous frequency identification for this little studied star, as well as to discuss the detected 
frequencies in the context of magnetic models and analyze the influence of the magnetic field on the pulsation. 
Methods. Using traditional Fourier analysis techniques on three independent data reductions, intrinsic frequencies for the star are 
identified. Theoretical non-adiabatic axisymmetric modes influenced by a magnetic field having polar field strengths Bp = 0-5 kG 
were computed to compare the observations to theory. 

Results. The high-precision data allow us to identify three definite intrinsic pulsation frequencies and two other candidate frequencies 
with low S/N. Considering the observed spacings, only one (50.95 /jHz) is consistent with the main sequence nature of roAp stars. The 
comparison with theoretical models yields a best fit for a 1.95 Mq model having solar metallicity, suppressed envelope convection, 
and homogenous helium abundance. Furthermore, our analysis confirms the suspected slow rotation of the star and sets new lower 
limits to the rotation period (P rot > 1 month) and inclination (;' > 30 + 10°). 

Conclusions. The observed frequency spectrum is not rich enough to unambiguously identify a model. On the other hand, the models 
hardly represent roAp stars in detail due to the approximations needed to describe the interactions of the magnetic field with stel- 
lar structure and pulsation. Consequently, errors in the model frequencies needed for the fitting procedure can only be estimated. 
Nevertheless, it is encouraging that models which suppress convection and include solar metallicity, in agreement with current con- 
cepts of roAp stars, fit the observations best. 
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1. Introduction 

Rapidly oscillating Ap (roAp) stars are classified as chemi- 
cally peculiar main sequence stars pulsating with periods rang - 
ing from 6-21 mins. Since the first discovery by iKurtzl (Il978h . 
more than 30 members have been identified and their pulsations 
are generally believed to be high-order low-degree magneto- 
acoustic modes. The roAp stars are characterized by strong 
global magnetic fields with polar strengths up to several kG. 

The discovery of a synchronous modulation of the pul- 
sation amplitude with the magnetic field strength led to the 
oblique pulsator model, which assumes the pulsation axis is 
aligned with the magnetic axis which is itself oblique to the 
axis of rotation. This concept naturally explains the observed 
frequency triplets as axisymmetri c dipolar modes modulated by 
the rotation period (Kurtzj 1982). This model was challenged 



by iBigot & Dziembowskil (|2002|) . who tookg into account ro- 
tational effects such as the centrifugal force and argued in fa- 
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vor of non-axisymmetric modes. As suming an axisymmetric 
mode aligned with the magnetic axis, ISaio & Gautschyl d2004t) 
found that the amplitude of the pulsation is confined to the mag- 
netic a xis, which was confirmed observationally by lKochukhovl 
d2006l) analyzing spectroscopic line profile variations of the 
rapidly rotating roAp star HR3831. 

The mechanism driving the pulsation in roAp stars is widely 
believed to be the ^-mechanism acting in the H ionization zone. 
Attempts have been made to explain the lack of low-order pul- 
sation modes despite the vicinity to classical <5Scuti pulsators 
in the Hertzsprung-Russell diagram (HRD) by both indirect and 
direct effects of the magnetic field. With the inclusion of grav- 
itational settling of helium in the outer layers of the star due 
to inhibited meridional circulation as well as suppression of en- 
velope convection at the poles due to the magnetic field, it be- 
came possible to driv e the high-order pulsations as observed 
dBalmforth et alj200ll) : however, some lower order modes were 
also excited in these models. 



I Said d2005l) concentrated on the direct effect of the magnetic 
field, showing that the decoupling of the pulsation into an acous- 
tic wave and a magnetic slow wave (which dissipates in the inte- 
rior) is sufficient to damp out all the low-order modes, while high 
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order modes remain unstable. Beside the excitation analysis, the 
direct effect on the observed frequencies has bee n intensively 
studied (e.g. JSaio & Gautschy||2004tlCunhall2006l) . 

It has been shown that magnetic coupling in the layers where 
the magnetic pressure becomes comparable to the gas pressure 
can significantly alter the pulsation frequency. This, in turn, 
causes non-negligible deviations from the classical asymptotic 
theory (lTassoullll980h . which predicts equally spaced frequen- 
cies of consecutive radial overtones and the same degree and 
is widely used to interpret high-order p-modes for, e.g., solar- 
type oscillations. It was shown that the observed deviations from 
these spacings depend essentially on the strength of the magnetic 
field as well as the field geometry, making accurate frequency 
determination a valuable tool in mode identification. 

The rapid pulsation of 10 Aql (HP 176232, HR 7167, V = 
5.91) was first discovered bv lHeller & Kramerl(ll990l) using high 
speed photometry with the 1-m telescope of the Mount Laguna 
observatory. From about 30 hrs of observations, they were able 
to identify three pulsation frequencies at /i,hk = 1.4360mHz, 
Ahk = 1.3854mHz and /^hk = 1.2393mHz all with am- 
plitudes < 0.5 mmag. Due to severe aliasing in their data, and 
the low signal amplitudes, they noted that the identifications are 
not unambiguous; however, they estimated from the observed 
spacings that the large frequency separation (hereafter sim- 
ply Ay ) could be 50.6^Hz. Shortly thereafter. iBelmonte et al.l 
(1 199 lb published IR photom etry which confirmed the detection 
bv iHeller & Kramerl (1 19901 ) with the signals in their amplitude 
spectra being close to the observed frequencies. Surprisingly, 
their observed amplitudes (obtained in the / filter) were in some 
cases about a factor 2-3 greater than in the Johnson B filter, and 
their high est peak was coincident with a 1 d alias of /3,hk- One 
year later. IBelmonte et ail (1 19921) returned to 10 Aql, but failed 
to redetect the signal by the same method. 

IR photometry in the early 1990's constituted the last pub- 
lished photom etric campaign of 10 Aq l. Focusing on abun- 
dance analysis, iRvabchikova et all d2000l) obtained high resolu- 
tion spectra of 10 Aql, estimating r e fj, radius and \ogg through 
synthetic line profile fitting. As for many other roAp stars, they 
detected a strong overabundance of the doubly ionized rare 
earth elements (REE's) Pr III and Nd III. Shortly thereafter, 
IKochukhov et a l. (2002) succeeded in detecting the pulsation 
from radial velocity variations of these lines while also confirm- 
ing the unusually low amplitude of the signal with peak ampli- 
tudes of 130-150ms~'. They also provided the first measure- 
ment of the mean mag netic field modulus (B) =1 .5 + 0.1 kG 
from Zeeman splitting. lHatzes & Mkrtichianl (120051) confirmed 
the pulsation in the radial velocities for several lines. They also 
detected highly variable amplitudes from line to line, but also 
from night to night for the same line, which was interpreted as 
either rotational modulation, growth or decay of modes on short 
timescales, or close multiperiodicity. As the MOST observations 
will show, the latter explanation was most likely the reason. 

In an effort to suggest 10 Aql as a potential target for th e 
space telescop e COR OT dBaglin & The COROT Team! 1 19981) . 
iBigot & Weissl d2002l) published the first and, so far, last at- 
tempt to model the oscillation spectrum of this star. Assuming 
an effective temperature of 8000 + 100K (based on photometric 
calibrations) and a luminosity of 21.4 + 0.6 Lq (obtained with 
the Hipparcos parallax), t hey found an agreement with the three 
frequencies published bv lHeller & Kramerl d 19901) and the large 
frequency separation expected for a magnetic field strength of 
0.8 kG. 

The rotation period of 10 Aql is unknown; however, it is 
believed to be very long based on magnetic field measure- 



Table 1. Basic data on the MOST observations of 10 Aql. 



Parameter 


Value 


Obs.Time [HJD-2451545] 


2371.08-2402.24 


Total timespan [d] 


31.16 


Duty cycle [%] 


98 


Exposure time [sec] 


10 


Sampling time [sec] 


20 



ments by Babcockl dl958l) . w hich showed a reversal of polar- 
ity in a timespan of 4 years. iPrestonl d 19701) measured a con- 
stant mean lon gitudinal field strength of ~ 500 G over 10 con- 
secutive nights TRvabchikova et alj d2005l) obtained 5 additional 
values over a timespan of one year, confirming the value mea- 
sured by Preston some 30 years earlier. Considering the mea- 
sured vsirn = 2.0 + 0.5kms _1 dKochukhov et al.ll2002l) as an 
upper limit, they concluded that the rotational period might pos- 
sibly be as long as hundreds of years. Since this, however, would 
contradict Babcock's observation and a pole-on configuration 
could not yet be excluded, there still was no securely measured 
lower limit for the rotation period of 10 Aql. 

All spectroscopic campaigns conducted on 10 Aql were 
based on observations spanning a couple of nights and were 
therefore unsuitable for confirming and precisely identifying the 
pulsation modes of the star. Although recent observations have 
shown the power of high resolution spectroscopy in probing the 
interior and eve n mapping the sur face structure of pulsation in 
roAp stars (e.g.. lKochukhovl2004l) . it is clear that a secure iden- 
tification of all pulsation modes is indispensable for a complete 
understanding of the effects observed (in both photometry and 
spectroscopy). Due to the high precision and long uninterrupted 
time coverage, MOST is uniquely suited for such a campaign 
to perfectly complement recent accomplishments in the field of 
roAp stars. 



2. Photometry and Data Reduction 

The MOS T (Microvariability and Oscillations of STars) space 
telescope (Wa lker et al.ll2003l) was launched in June 2003 to per- 
form high precision uninterrupted photometry of bright stars. It 
is equipped with a Maksutov telescope with a 15 cm aperture, a 
custom broadband filter (350 - 700 nm) and, in its original de- 
sign, two CCD's (one for science, the other for startracking). 
Due to its sun-synchronous, polar orbit with an orbital period 
of ~ 101 mins, MOST is able to measure stars for up to two 
months without interruption, making it uniquely suited to ob- 
tain high precision frequency determinations free from aliasing 
effects. Early in 2006 the tracking CCD system failed due to a 
particle hit. Since then, both science and tracking have been car- 
ried out with the Science CCD system. 

The star 10 Aql was observed by MOST from June 29 - July 
30 2006 as a primary science target. Because of its moderate 
brightness (V = 5.91) it was observed in th e open field (corre- 
sponding to traditional CCD photometry, see IWalker et al.ll2003l 
for details). A single exposure consisted of 20 added 0.5 sec im- 
ages for a total exposure time of 10 sec at a sampling interval of 
20 sec. 119724 exposures were acquired with two interruptions 
due to difficulties with the satellite electronics, to yield a duty 
cycle of 98%. 

Due to parasitic stray light caused mostly by Earthshine 
scattered into MOST's focal plane, special data reduction tech- 
niques were developed to remove instrumental artifacts present 
in the raw data. In the case of 10 Aql, three different indepen- 
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dent reduction techniques were applied. A newly developed re- 
duction software develo ped for MOST open field photometry 
dHuber & Reegenl l2008) is b ased on a techniqu e originally de- 
veloped for Fabry imaging dReegen et al.l 1200 6). in which cor- 
relations of intensities of pixels illuminated by starlight and sky 
intensities are calculated over time. Performing a correction of 
a linear regression to the intensity-intensity diagrams, stray light 
effects are effectively eliminated from the data. To be able to 
apply this technique to open field photometry, the reduction ap- 
plies an image shifting routine which aligns all images to com- 
mon pixel coordinates using bilinear interpolation of sub-pixel 
intensities, allowing a definition of a fixed aperture over time in- 
dependent of minor satellite jitter effects. The software allows 
a decorrelation of the mean as well as pixel intensities with the 
option of using higher order fits (depending on the complexity of 
the correlation). It additionally implements an identification and 
correction of pixels affected by cosmic ray hits as well as elimi- 
nation of bad quality fr ames caused by pointing instabilities (see 
iHuber & Reegen1l2008l for details). 

An independent reduction using the ope n-field subraster wa s 
performed by JR and is described in detail bv lRowe et al.l(l2006l) . 
This technique additionally implements flatfield and dark current 
corrections. Since there are no moveable parts on the satellite to 
obtain such calibration frames in the traditional sense, flatfield 
correction is performed by making use of the fact that during 
high stray light phases individual CCD subrasters are almost uni- 
formly illuminated by scattered Earthshine. The software uses a 
combination of PSF fitting and aperture photometry to obtain the 
count rates of the star. Stray light influences are corrected similar 
to the reduction described above from a fit to mean star and sky 
intensities. 

With the loss of the MOST startracking CCD another reduc- 
tion capability was implemented. MOST acquires count rates for 
several guide stars, which are reduced on-board by simply sub- 
tracting a mean sky value after applying a given threshold. Since 
guide stars and science targets are on the same CCD, the guide 
star values are also returned for all science targets (in the case of 
10 Aql, every 0.5 sec). The extracted on-board processed inten- 
sities were additionally reduced by MH by applying the decor- 
relation technique with other (constant) guide stars. While the 
accuracy of the count rate determination is certainly not of the 
same quality as in the actual subraster images, the big advantage 
of this reduction is that intensities are on-board reduced after 
each sub-exposure, avoiding possible smearing effects of stray 
light effects in stacked images. 

All three reductions were performed independently and 
yielded different effective duty cycles (ranging from 63 - 95 %) 
depending on the outlier rejection philosophy used. Each reduc- 
tion showed several strengths and weaknesses, and all three were 
used to secure the results obtained from the data. FigureQ]shows 
the reduced light curve of 10 Aql, in which all signals with fre- 
quencies < 5 d" 1 corresponding to instrumental trends have been 
removed. The data show an increased scatter towards the end of 
the run. Nevertheless, the bottom panel clearly shows the detec- 
tion of roAp pulsation. 



3. Frequency Analysis 

3.1. Frequency Identification 

Frequency anal ysis was performed usin g the programs SigSpe c 
dReegenl |2007|) as well as Period04 dLenz & Bregerl |2004|) . 
SigSpec computes a quantity called spectral significance based 
on a false-alarm probability by comparing the signal to white 
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Fig. 1. Reduced light curve of 10 Aql. The bottom panels show a 
zoom of the light curve, at approximately HJD-245 1545 2396.9- 
2397.0 with a running average overplotted. 



noise, and incorporating frequency, amplitude and phase infor- 
mation in the calculation. The program has proved to be ex- 
tremely powerful for a statistically clean evaluation of the sig- 
nal in Fourier space, and especially in the case of the high fre- 
quencies for the roAp stars where colored noise (intrinsic or in- 
strumental) can be assumed to be negligible, the spectral signif- 
icance gives important additional information compared to the 
more traditional S/N estimate. 

Figure [2] (top panel) shows a Fourier amplitude spectrum 
of the reduced 10 Aql data in the frequency region of inter- 
est. With the dashed lines indicating the pu blished frequency 
values of the two largest amplitude peaks by iHeller & Kramerl 
dl990h . the MOST data clearly show that indeed their frequen- 
cies were misidentified as 1 d -1 (= 0.0116mHz) aliases. In fact, 
for both cases the correct frequency lies 1 d -1 higher. Moreover, 
the MOST data enables us to resolve an additional, as of yet 
undetected, frequency around 1.427 mHz at high S/N. Thanks 
to the clear spectral window showing only the aliasing effects 
spaced at MOST's orbital frequency (14. 19 d -1 = 0.16 mHz), the 
identification of these three frequencies is unambiguous. They 
appear consistently in all three reductions with comparable sig- 
nificances. To ensure the intrinsic nature of this signal, the mid- 
dle panel of Figure[2]shows the amplitude spectrum of the bright- 
est star simultaneously observed (HD 175922, V = 7.21) and ad- 
ditional peaks that were identified using the same approach with 
SigSpec to the light curve constructed with background readings 
of the 10 Aql subraster. Clearly, no instrumental power can be 
identified in the vicinity of f\ , ft or . We therefore definitely 
consider these frequencies to be the result of intrinsic stellar pul- 
sation. 

Apart from the three obvious peaks, 10 Aql does not reveal 
any eye-catching power above the noise level of about lOppm. 
We do detect two peaks in the spectrum which were identified by 
SigSpec with a spectral significance of 6, corresponding roughly 
to a S/N > 4. These frequencies are marked in Figure [2] by the 
labels fi and /j. Unfortunately, these peaks are only identified 
in the reduction showing the lowest noise level in this region. 
The comparison data (middle panel, Figure |2]i does not reveal 
any obvious power in this frequency region; however, one must 
note that the amplitudes of f\ and /j are at about the noise level 
of that in the fainter comparison star. Clearly, we are not able to 
infer from the data available if these peaks are indeed intrinsic 
or not. Nevertheless, we consider them as candidate frequencies 
and included them in the subsequent analysis. 
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Fig. 2. (Top panel) Amplitude spectrum of 10 Aql. Solid lines indicate the position of published frequencies by iHeller & Kramerl 
dl990j) . Dotted lines show the position of 1 d _1 sidelobes of harmonics of MOST's orbital frequency (i.e., instrumental artifacts). 
The insert displays the spectral window. (Middle panel) Inverted amplitude spectrum of the brightest guide star simultaneously 
observed (HD 175922, V = 7.21), shown at the same scale. Dashed-dotted lines indicate the position of frequencies identified in the 
background readings of the 10 Aql subraster. (Bottom panel) Frequency region where Ahk was reported after prewhitening of f\, 
f2 and /3 (the position again indicated by a solid line). 



Finally, the bottom panel of Figure [2] shows the fre- 
qu ency region where an a dditi onal mode has been reported 
by IHeller & Kramerl (1 1990b and lBelmonte et al.l d!991l) . We do 
not detect any significant power at A hk or an V surroundin g 
1 d -1 aliases thereof (as suggested by iBelmonte et al.l 119911) . 
Considering that Ahk was the lowe st amplitude frequency in 
the dataset of IHeller & Kramerl (1 19901) . one might conclude from 
this result that their detection was indeed an instrum ental arti- 
fact. Nevertheless, the fact that lBelmonte et ail d 19911) also found 
a signal at almost exactly the same value (in their case, the high- 
est amplitude peak) makes this explanation more unlikely. On 
the other hand, their follow-up observations one year later did 
not show evidence for any detec table signal at all using the same 
methods (IBelmonte et al.fll992l) . 

If the detections were real, we could conclude that the lack 
of signal in the MOST data is in fact evidence for a limited mode 
lifetime on the order of some decades, possibly caused by damp- 
ing of the mode by intrinsic changes in the star, or a different 
orientation of the star at the time of MOST observations due 
to rotation over the past 15 years. Due to the uncertainty of both 
earlier detections, however, these conclusions must remain in the 
realm of speculation at this point. 



Apart from the frequency region just discussed, the 10 Aql 
data do not reveal any variability at lower (or higher) frequen- 
cies with a S/N > 4 which cannot be attributed to instrumen- 
tal artifacts. This is consistent with previously published results 
demonstrating 10 Aql to be a slow rotator (i.e., no frequency 
multiplets due to rotational modulation are observed) and the- 
oretical predictions that no low-order p-modes are excited in 
roAp stars. We are, furthermore, not able to detect harmonics of 
any definite frequencies f\, A or fa, contrar y to what has been 
found for the MOST observations of y Equ dGruberbauer et all 
120081) . Table|2]lists all identified intrinsic and candidate frequen- 
cies in the MOST 10 Aql dataset. Note that the periods of the 
stated frequencies correspond to 11.3 - 12.2 minutes and that 
errors of the frequency, amplitude and phase values have been 
computed incorporatin g the spectral significance as described in 



computed incorporating 
iKallinger etai1(l2008l) . 



3.2. Amplitude Modulation 

Apart from mean light variations which are believed to be due 
to abundance spots on the surface of Ap stars, amplitude mod- 
ulation of an intrinsic frequency might present clues about the 
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Table 2. List of frequencies identified in the MOST 10 Aql ob- 
servations. Candidate frequencies are denoted /; and /j. The 
phase values shown are calculated as sin(27r(/f + 6)) with the 
zero point corresponding to JD 2451545. 



Id 


/■[mHz] 


A[mmag] 


<9[rad] 


sig 


fi 


1.44786 ± 0.00003 


0.17 ± 0.01 


0.2 + 0.3 


135 


h 


1.39691 ± 0.00003 


0.15 ± 0.01 


0.8 + 0.3 


101 


h 


1.42709 ± 0.00004 


0.12 ± 0.01 


0.7 ± 0.4 


68 


fi 


1.3662 ± 0.0001 


0.03 ± 0.02 


0.3 ± 1.2 


6 


fi 


1.4686 ± 0.0001 


0.03 ± 0.02 


0.9 ± 1.2 


6 



60 



40- 




-40 - 



5 10 15 20 25 

Time [d] 

Fig. 3. Amplitude modulation of f\ as determined by time re- 
solved frequency analysis. The different lines correspond to the 
three reductions performed for the dataset (solid: Direct Imaging 
(Vienna); dashed: Direct Imaging (UBC); dashed-dotted: Guide 
Star), and each line shows a smoothed curve of the individual 
values. The ordinate is the relative amplitude change. The error 
bar indicates the mean error for all three reductions for a single 
datapoint. 



rotation period. Following the oblique pulsator model, the mag- 
netic poles (and hence the pulsation poles) would wander in and 
out of sight as the star rotates, modulating the pulsation ampli- 
tude with exactly that period. As the 31.2 days of MOST obser- 
vations constitute the longest continuous dataset ever obtained 
for 10 Aql, we are able to set a new limit to the so far unknown 
rotation period of this star. 

The modulation was investigated by time resolved frequency 
analysis using SigSpec Subsets were chosen long enough to re- 
solve the frequency component with a secure detection limit 
(in this case 3.2 d), with stepwidths equal to one pulsation pe- 
riod. This exercise was performed for all three intrinsic frequen- 
cies and all three independent reductions. The results for f\ are 
shown in Figure[3] While the amplitude does change up to 25 % 
throughout the run, none of the variations are consistent within 
the different reductions. Similar results have been obtained for 
fz and /j. We conclude that the variations seen are purely due 
to different data and reduction qualities at different observation 
times, and that no significant amplitude modulation is present in 
the 10 Aql data. Considering the timespan of 31.2 days of ob- 
servations, we can therefore rule out that the star rotates with 
a period shorter than about a month (assuming that the star is 
not seen pole-on, in which case no modulation would be de- 
tected even for a rapidly rotating star). This conclusion is also 
confirmed by time-resolved spectroscopy which was obtained si- 
multaneously to the MOST photometry spanning over 20 days, 
yieldi ng constant line intensities within 0.5 % dSachkov et alj 
120081) . Both results are consistent with earlier observations that 



10 Aql is a slow rotator and set a new lower limit for its rotation 
period. 

4. Asteroseismic Analysis 

4. 1 . Frequency Spacings 

A fundamental observational parameter in stars pulsating with 
high-order acoustic modes is the large frequency separation Av, 
describing the regular spacing of consecutive radial overtones n 
with same degree t. Since it is proportional to the square root of 
the mean density of the star, it allows an estimation of astrophys- 
ical parameters such as luminosity and effective temperature. 

Obviously, the intrinsic frequencies observed in 10 Aql do 
not reveal a regular spacing. The observed spacings are as fol- 
lows: 

/i - f 2 = Ai = 50.95 ± 0.03 /iHz 
h ~h = A 2 = 30.19 + 0.03 /iHz 
/i - f 3 = A 3 = 20.76 ± 0.03 yuHz 

Considering roAp theory as mentioned in the introduction, 
this result is not surprising, since it has been shown that the mag- 
netic field can significantly perturb the observed frequencies and 
therefore the spacing. On the other hand, it might be possible 
that what we observe are in fact modes of different degree I. 
It is, however, interesting to note that f% — f\ ~ A2 as well as 
fj-fi ~ A3 . This might give rise to the assumption that either f\ 
or fi are perturbed by the magnetic field, which, if those values 
were slightly shifted, give four consecutive frequencies spaced 
by either A2 or A3 . 

To test these speculations, we used the published HRD 
(Hertzsprung Rus sell Diagram) position by | Matthews et all 
d 19991) as well as iKochukhov & Bagnulol (120061) . Thev calcu- 
late the luminosity based on the measured Hipparcos parallax 
and derive an effective temperature based on photometric cal- 
ibrations. The frequency spacing, Av, is directly connected to 
physical properties which determ ine the position of the star in 
the HRD, as given for example by [Matthews et al.ld 19991) via 

Av = (6.64 + 0.36) x 10" 16 M°' 5 7; e 3 ff L" ' 75 Hz , (1) 

where M is given in solar masses, L in solar luminosities and 
r e ff in Kelvin. Hence, it is possible to identify models show- 
ing a constant Av and compare them with the previously deter- 
mined values. The result is shown in Figure|4] In this theoretical 
HRD, evolutionary tra cks using the Yale Stellar Evolution Code 
dGuenther et al.l 1 992h are plotted, together with lines of constant 
Av of 50 and 30 /iHz, respectively. 

Clearly, any large frequency separation considerably less 
than 50 /iHz would contradict (for the given temperature and lu- 
minosity region) the current understanding of roAp stars as ob- 
jects between ZAMS and TAMS. Note that although these evo- 
lutionary tracks do not take into account any effects of the mag- 
netic field, frequency jumps on the order of several tenths of //Hz 
are expected to occur only at discr ete frequency values rather 
than constantly for each frequency dCunhall2006l) . More likely, 
consecutive shifts are on the order of a few /iHz. It is therefore 
not possible that the magnetic field perturbations mimic a Av 
which is different from that expected from the physical parame- 
ters of the star. In other words, equation dTJ is still a valid approx- 
imation even for modes perturbed by a magnetic field, except for 
large jumps at certain overtones. It is interesting to note that this 
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lo S T e rr 

Fig. 4. Theoretical HRD showing ev olutionary tracks wit h 
(X,Z) = (0.71,0.19) and a = 1.8 dGuenfher et al.1 [l992b . 
Numbers next to the tracks show the mass in solar units. The 
dotted lines indicate contours of constant Av. Previously deter- 
mined parameters for 10 Aql are shown. 

Table 3. Summary of models used for the analysis. 



Id 


M/Mq 


(X,Z) 


He depl 


Conv 


S190 


1.90 


(0.700,0.020) 


yes 


no 


S195 


1.95 


(0.700,0.020) 


yes 


no 


S200 


2.00 


(0.700,0.020) 


yes 


no 


S205 


2.05 


(0.700,0.020) 


yes 


no 


H195 


1.95 


(0.700,0.020) 


no 


no 


H200 


2.00 


(0.700,0.020) 


no 


no 


H205 


2.05 


(0.700,0.020) 


no 


no 


M200 


2.00 


(0.695,0.025) 


yes 


no 


C195 


1.95 


(0.700,0.020) 


no 


yes 



value for Av was act ually already suggested a fter the first photo- 
metric campaign bv lHeller & Kramerj(ll990l) . Although less ac- 
curate and based on wrong frequency identifications, their first 
asteroseismological interpretation of 10 Aql might turn out to be 
correct after all. 

4.2. Model Frequencies 

To further quantify the observations of 10 Aql, we performed a 
comparison with models which were calculated from the zero- 
age main sequence to the end of the main sequence for masses 
ranging from 1 .90-2.05 Mp) in steps o f 0.05 Mq using OPAL 
opacity tables (llglesias & Rog ers 1996). In standard models, an 
initial chemical composition (X, Z) = (0.70, 0.02) is adopted and 
envelope convection is suppressed. Helium abundance is taken 
as Y = 0.01 + 0.27(^2 + X3), where Y is the helium mass fraction 
and X2 and represent fractions of singly and doubly ionized 
helium, respectively. For standard models, helium is depleted 
at the surface (Y$ = 0.01) due to gravitational settling in the 
outermost layers of the star. Additionally, some non-standard 
models were computed which included metal-rich composition 
(Z = 0.025), envelope convection or assumed no helium deple- 
tion near the surface (Fs = 0.28). Table [3] lists all the stellar 
models used in our analysis. 

Theoretical non-adiabatic frequencies of axisymmetric (m = 
0) high-order p-modes under the influence of a dipole magnetic 



field were calculated using the method described in lSaiol d2005l) . 
Outer boundary conditions are imposed at an optical depth of 
10~ 3 . A reflective mechanical condition {Sp/p — > constant) was 
used although for models with Y$ = 0.01, the acoustic cut-off 
frequency is close to the frequencies observed for 10 Aql (for 
models without helium depletion the critical frequency is higher 
than the observed frequencies, therefore justifying the use of a 
reflective boundary condition). Due to the neglect of possible 
kinetic energy leakage at the outer boundary, and the use of a 
diffusion approximation for radiative energy transport, the pre- 
diction for the stability of the modes in these models may not be 
very reliable. 

Eigenfunctions of theoretical frequencies are expanded in 
terms proportional to spherical harmonics Y"' =0 (9, (p). The char- 
acter of the latitudinal amplitude variation is designated by £■&, 
the I value of the expansion component having the largest ki- 
netic energy. The value of 1% can change as the eigenfunction of 
a mode is modified by the effects of the magnetic field. In fact, a 
£b = mode sometimes changes to a -2 mode or vice versa. 
Generally, frequencies of — 0, 3 modes are strongly affected 
by a magnetic field, while i% = 1,2 modes tend to be less af- 
fected. The latter modes can therefore be used to determine the 
position in the HR diagram and the former modes to estimate the 
strength of the magnetic field. 

Theoretical frequencies were calculated for the HRD posi- 
tion estimated by previous observations as well as the observed 
frequency spacings and for magnetic field influences ranging 
from B P = 0-5 kG. 

Because of the necessary simplifications, the models are 
far from representative of roAp stars in detail. For example, 
the magnetic field is assumed to be a simple centered dipole. 
Consequently, errors of model frequencies needed for the com- 
parison with observational results can only be estimated. An ad- 
ditional severe restriction in the models is that only axisymmet- 
ric modes are considered. In the case of 10 Aql, however, this 
assumption can be justified by the fact that it has been repeat- 
edly shown that this star has a very long rotation period possibly 
in the order of years, in which case rotationally perturbed modes 
are negligible. As has been shown in section [3~2l this is also sup- 
ported by MOST photometry. 

4.3. Model Fitting 

Due to the fact that the calculation of model frequencies is 
very time consuming, the resolution of the original grid of 
models was enhanced b y linear interpolation as described by 
iGruberbauer et al.l ((2008) in the application to the MOST obser- 
vations of y Equ. The stepsize was decreased to 0.02 kG in Bp as 
well as 0.00001 in logL/Lg and log T e g. 

Observed frequencies were compared to theoretical ones 
using the x 1 statistics test as described by iGuenther & Brownl 
(120041) and defined as 

2 1 l, v obs,i _ V mo d,iJ ... 

X = TT ), — 2 ' (2) 

/— 1 obs.t mod.t 

where v b s ,i and v mo dei,i are me observed and corresponding 
model frequency, and cr obs i and cr mod i are the uncertainties of 
these values, respectively. N denotes the total number of fre- 
quencies fitted. For the uncertainty of model frequencies a value 
of 0.2 //Hz was assumed. Using this formalism, a set of observed 
frequencies is considered to match the theoretical frequencies 
within the errors for values of ^ 2 < 1 . We note that with this 
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H195(f/f 2 /f 3 ) 



H195(f/f 2 /f 3 /f.) 




49 50 51 52 53 1 2 3 4 5 
Av [uHz] B [kG] 




0.01 0.02 0.03 0.04 
Frequency mod 0.04974 



Hi95(f/f 2 /f 3 /f.) 



49 50 51 52 53 1 2 3 4 5 
Av [nHz] B [kG] 



50 51 52 53 1 2 3 4 5 
Av [uHz] B [kG] 




0.01 0.02 0.03 0.04 0.05 
Frequency mod 0.04979 



0.01 0.02 0.03 0.04 
Frequency mod 0.05049 [mHz] 



Fig. 5. Graphical presentation of the model fitting results for different frequency combinations in the H195 model. (Top panels)^ 2 
values as a function of the parameter space covered by the models in Av (i.e., mass, effective temperature and luminosity) and Bp, 
respectively. (Bottom panels) Echelle diagrams of models with the lowest^- 2 . Spherical degrees ^ — (triangles), — 1 (squares), 
{b =2 (circles) and £r = 3 (diamonds) are displayed, with observed MOST frequen c ies lab eled and marked by plus signs. The cross 
in each plot shows the position of a frequency value published by iBeimonte et al.l dl99ll) at 1.247 mHz. Both errors for observed 
frequencies and theoretical frequencies are slightly smaller than the symbol sizes. 



formulation we do not draw any conclusions about over- or un- 
derestimation of errors (e.g., a reduced ^- 2 <k 1 value would, 
in the strictly statistical sense, indicate an overestimation of er- 
rors), but rather solely use the relative x 2 values to differentiate 
between the quality of the fits. 

Since the intrinsic nature of the frequencies fi and /j cannot 
be confirmed beyond doubt, model fitting was performed with all 
combinations of definite and candidate frequencies. The results 
are summarized in Table |4] As can be seen from the^- 2 values, 
the observed frequencies fit the theoretical ones well, if not to 
say too well considering the rather crude magnetic stellar mod- 
els. Interestingly, when using all five frequencies, no fits within 
the errors can be found. This naturally leads to the conclusion 
that, assuming the theoretical frequencies are correct, either fi 
or fj must not be intrinsic to the star. On the other hand, if only 
the definite frequencies are used, all models fit within the er- 
rors, showing that an additional candidate frequency is needed 
to further constrain the models. Therefore, the focus of attention 
was brought to the combinations including only one of the two 
candidate frequencies. 

Although the diversity of x 2 values for all models is quite 
weak, one can see that (apart from the combination using only 
definite frequencies), some non-standard models seem to clearly 
yield worse fits. This is true for models including envelope con- 
vection as well as models with higher metallicity Z. The former 
result is consistent with roAp theory, showing that suppressed 
envelope convection is a necessary condition for driving pulsa- 
tion in these stars (e. g„ iBalmforth et alJ[200U ISaioll2005l) as has 



Table 4. Model fitting results for each combination of definite 
and candidate frequency set. Bold numbers show the best fitting 
results for a combination. 



Model 






x 1 




id 


(/1//2//3) 




(fi/fz/fi/fi) 


<Jv'f:>f ■•>/„>/■) 


S190 


0.27 


1.51 


0.94 


6.35 


S195 


0.13 


0.55 


1.33 


7.96 


S200 


0.10 


0.38 


0.69 


2.10 


S205 


0.26 


0.46 


2.63 


2.61 


H195 


0.05 


0.13 


0.29 


4.54 


H200 


0.02 


2.62 


2.05 


3.50 


H205 


0.13 


1.40 


4.71 


4.08 


C195 


0.11 


0.84 


2.33 


8.32 


M200 


0.30 


5.20 


4.33 


15.2 



also been found by lGruberbauer et al l (120081) for the MOST ob- 
servations of y Equ. 

Considering only frequency combinations where fits with 
X 1 < 1 can be found, it can be seen that the best fits consistently 
occur for models without helium depletion. Examining the best 
fit for the set of definite frequencies (H200), it was found that the 
corresponding polar magnetic field value lies at 0.7 kG. While 
thi s value seems t o be in ag reement with the values derived 
by lBigot& Wei ss (2002) and iKochukhov et alj d2002l) . it must 
be noted that these authors considered the mean magnetic field 
strength, whereas the pulsation models test the polar magnetic 
field, which is expected to be considerably higher. Therefore, 
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this result woul d stand in direct contrad iction to the observed 
value inferred bv lKochukhov et all d2002l) . Furthermore, the two 
primary frequencies are identified as radial ({% = 0) modes in 
this model, which contradicts evidence found hitherto for roAp 
stars as being non-radial pulsators at least for large amplitude 
frequencies. For these reasons, we conclude that also for the 
combination of /1//2//3 the HI 95 model, yielding a higher^ 2 
value (but still < 1) and a polar magnetic field value of ~ 3.7 kG, 
presents a more reasonable solution. 

In the top panels of Figure[5] x 2 1S shown as a function of Av 
and Bp for the HI 95 model. To present the results more easily, 
r e ff , L and M have been combined as Av according to Equation 
Q] Including fi, the regions showing low^- 2 reduce dramatically 
(as expected), with only one region in Av and Bp having x 2 < 
1. The same is valid for the combination including fi, which, 
however, does not restrict the parameter space as much as the 
former combination, with the best fit lying at a slightly higher 
Bp and higher Ay. 

As shown bv lSaio & Gautschvl d2004 and lSaiol d2005[) . there 
are certain values of Bp where the damping of acoustic waves 
due to magnetic slow waves reaches its maximum. In some cases 
the expansion of eigenfunctions used in the method does not 
converge, resulting in gaps in the theoretical frequencies. Note 
that the region between Bp = 2.3-3 kG shows such a gap for our 
set of frequencies. 

In addition to t he frequencies observed b y M OST, we test the 
values provided bv lHeller & Kramerl d 19901) and lBelmonte et all 
d!991l) for the possible (but undetected in MOST data) additional 
frequenc y around 1.24mHz. Ind eed, one of the three values pub- 
lished bv lBelmonte et al.l dl99ll) (1 .247 mHz) would fit the same 
sequence of fi and f2, as shown in Figure[5]by the cross in the 
corresponding Echelle diagram. Note that this is not the case 
for the combinations with fi and that for the combination us- 
ing fi, fewer theoretical modes are present since the fit lies very 
close to the critical Bp introducing a gap in the parameter space. 
Although this result is encouraging, we do not feel comfortable 
basing a model fit on a frequency that is not observed by MOST. 
Nevertheless, this is another indication that 10 Aql indeed might 
have pulsated with an additional, slightly longer period in the 
early 1990's. 

It might appear questionable to prefer the quality of one fit 
over the other when both show a^ 2 < 1, i.e. both solutions fit 
within the errors. We argue that this is reasonable with the fol- 
lowing example: the calculation of the minimum mean devia- 
tion of observed to theoretical frequencies, which does not in- 
corporate any error values, gives (e.g., for the combination of 
certain frequencies) a value of 0.03 //Hz for the H195 model. 
The corresponding standard model S195 (also showing a^ 2 < 1 
close to the H195 model, see TableHJ) yields a mean deviation of 
0.07 //Hz, meaning that the expected values of both theoretical 
and observed models are closer together for the model showing 
the lower x 2 value. If we now assume that the error distribution 
is Gaussian, it is clear that fits with a lower mean deviation from 
the expected values imply a higher probability that the solution 
is in agreement with the real (unknown) values. 

Despite of this argument, it is clear from the results that there 
are not enough definite frequencies in the MOST 10 Aql dataset 
to obtain a unique solution. Also, it is not possible to clearly 
rule out whether fi and fi are artifacts or not. Indeed, tests have 
shown that if random values or frequencies detected in simulta- 
neously observed stars with a comparable significance to fi or fi 
are added to the set of definite frequencies, in about 50 % of the 
cases a fit showing a^ 2 < 1 can be found for the HI 95 model. 
This result demonstrates the fact that such a small quantity of 
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Fig. 6. Flux perturbation as a function of co-latitude (cos = 
1 referring to the poles) for the best fitting non-interpolated HI 95 
model (log r eff = 3.8883, logL/L©=1.2473). Note that f 2 (4 = 
2, n — 27) is not specifically shown since the variation coincides 
with fi . 



observed frequencies does not allow one to derive unambiguous 
results for this star and all conclusions based on these fits should 
be treated with caution. Nevertheless, it can be concluded from 
the attempts presented here that the most reasonable fits seem 
to occur consistently for one stellar model, and that the results 
including fi seem to constrain the parameter space better for this 
model. 

The bottom panel of Figure[5]shows the Echelle diagrams of 
the models with the lowest x 2 , together with the observed fre- 
quencies and a mode identification. Clearly, it is quite danger- 
ous to derive a mode identification considering the low x 2 val- 
ues given in Table[4]for other models. Nevertheless, we note that 
for the best solution including fi the mode identification {is/n) 
would be as follows: fi = (2, 28), f 2 = (2, 27), f 3 = (1, 28) and 
fi = (0, 26). If fi is included, the modes would be identified as 
fi = (1,28), / 2 = (l,27),/ 3 = (3, 27) and fi = (2,28). Figure © 
shows, for the closest non-interpolated model to our best fit, the 
flux perturbation as a function of the co-latitute at the phase of 
maximum light for each mode. It can be seen that for all modes 
maximum light is concentrated towards ~ 0.7, and that the 
radial mode exhibits larger amplitudes towards the equator. 

4.4. Comparison with previous results 

Although the results of fitting are not unambiguous, we can com- 
pare our best fitting model with previously determined parame- 
ters derived using non-asteroseismological observables. The re- 
sult of this comparison is shown in a theoretical HRD in Fi gure 
U\ In addition to the resu l ts deri ved by lMatthews e t al. ( 199^) and 
iKochukhov & Bagnulol (120061) . we plot the position of our best 
model fit using the combinations with fi (filled diamond; note 
that the combination with fi only slightly changes this position) 
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Table 5. Astrophysical parameters derived for 10 Aql based on 
calculations using the observed spacing and model fits (using 
fi). Bold face numbers show parameters based on estimations 
and non-asteroseismic observables, respectively. Errors for the 
model fit have been estimated by considering only grid points 
fits with x 1 < 1 ■ The model mass is a fixed parameter. 



as well as an estimate using the observed large frequency spacing 
in combination with non-asteroseismological parameters by ap- 
plying equation ([TJ. The most reliable effective temp erature es- 
timation for 10 Aql is from lRyabchikova et alJ (12000) who fitted 
line profiles to the spectral energy distribution, yielding T e s = 
7550 + 150 K. The mass was assumed to be 1.95 + 0.05 Mq. 
Using these combined parameters (partially) based on astero- 
seismological observables, we also derive R/Rq = 2.3 + 0.2 
and \ogg = 3.99 + 0.06 , both values being in ag reement with 
previous determinations dRvabchikova et aLll2000b . 

There are two notable issues in this comparison. First, 
one can see that there is a significant discrepancy between 
the luminosity derive d from asteroseismology and Hipparcos. 
iMatthews et all (1 19991) were the first to perform a comparison of 
the parameters derived by both methods for a sample of roAp 
stars, including 10 Aql (using the large spacing suggested by 
iHeller & Kramer! 1 19901) . In general, they found a good agree- 
ment between the two methods; however, they noted that sys- 
tematically most Hipparcos luminosities were lower than those 
ones derived from Av, except for 10 Aql. Arguing the other way 
around, they also compared effective temperatures derived us- 
ing Av and the parallax and compared those with temperatures 
derived from HyS photometry, concluding again a good agree- 
ment but "In only one case (10 Aql) does Hipparcos predict a 
higher temperature." If we assume that indeed the measured par- 
allax is inaccurate and 10 Aql should follow the general trend, 
then a lower effective temperature and hence a lower luminos- 
ity (in agreement with our results) would be needed. Although 
this explanation is cert ainly tempting and the deviation noted by 
IMatthews et al.l (1 19991) for 10 Aql is intriguing, the ambiguous 
nature of our derived fits makes us cautious about coming to any 
definite conclusions. 

Second, the comparison shows that the effective temperature 
derivation based on spectroscopy tends to cooler values, while 
photometric calibrations predict a high er temperature (it has to 
be not ed though that the investigation bv lKochukhov & Bagnulol 
J2006h was aimed solely to study Ap stars in a consistent manner 
rather than deriving an exact value for T e s). While the asteroseis- 
mic value lies in between, it can be seen that only a smaller por- 
tion of our models reach into the low temperatures as predicted 
by spectroscopy (see Figure[7]l, and that our best fitting solution 
does tend to cooler models, considering the grid investigated in 
this work. If the observed large frequency separation is indeed 
correct, such a low value for T e g would present another indi- 
cation for an overestimated luminosity based on the Hipparcos 
parallax. 
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Fig. 7. Theoretical HRD showing previously determined posi- 
tions of 10 Aql using calibrations based on photometry and the 
Hipparcos parallax (full circle and triangle) as well as the re- 
sults derived in this work. Evolutionary tracks correspond to the 
models presented in section 14.21 (the numbers next to them in- 
dicating the corresponding mass). Note that models excluding 
He depletion and including envelope convection show no signif- 
icant deviation from these tracks. The box shows the position 
derived assuming Av = 50.95 ^Hz, using effective temperatures 
derived from spectroscopy and an assumption of the mass to be 
1.95 + 0.05 Mq. Small crosses show the position of the mod- 
els listed in Table [3] with the filled diamond indicating the best 
model fit using a combination including f\. 



5. Conclusions 

MOST photometry confirms the rapid pulsation of 10 Aql, pro- 
viding for the first time an unquestionable identification of three 
pulsation modes in this star. Furthermore, we identify two can- 
didate frequencies with too low a S/N for a clear detection. We 
do not identify any signal at a frequency previously published in 
two independent ground-based datasets. We, nonetheless, find 
that this mode fits well with our present solution and conclude 
that it might be evidence for finite mode lifetimes in roAp stars, 
if the detections are real. Our results confirm previous findings 
that 10 Aql pulsates with very low amplitudes, and that if addi- 
tional modes are present their amplitudes must lie significantly 
below a level of about 40ppm. 

Based on the detected signal we can conclude that pre- 
viously detected amplitude modulation in spectroscopy by 
lHatzes & Mkrtichianl d2005l) was most likely due to beating of 
unresolved intrinsic frequencies. This is confirmed by simul- 
taneous spectroscopic monitoring, which shows that the ob- 
served amplitude mo dulation in 10 Aql is a pure beating effect 
(ISachkov et alj|2008l) . 

The MOST light curve presents the longest continuous 
dataset ever obtained for this star. A time-resolved analysis 
of the amplitudes for the detected frequencies yields no sig- 
nificant modulation. We conclude that in case the star is not 
seen pole-on, 10 Aql must rotate on times cales longer than a 
mont h, confirming previous findings (e.g., iRvabchikova et all 
2005 ). This is supporte d by the measured vsin/ = 2.0 ± 0.5 
dKochukhov et alj2002l) . which, together with the radius derived 
in this work, yields an approximate lower limit for the inclina- 
tion of; > 30+ 10°. 



Parameter using Av obs model fit J* 

M/Mq 1.95 ± 0.05 L95 g> ' 

r cff 7550 ± 150 7734.5 ± 0.5 3 

L/L 15.6 ± 1.7 17.6885 ± 0.0005 
Bp [kG] - 3.07 ± 0.03 

R/Rq 2.3 + 0.2 2.3488 + 0.0003 u 

logg 3.99 + 0.06 3.9860 + 0.0001 
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Our attempt to fit the observed frequencies to theoretical 
models does not yield an unambiguous solution. We note, how- 
ever, that models including envelope convection or higher metal- 
licity give a poorer fit. For all frequency combinations fitting 
within the errors, a 1.95 Mq model without helium depletion 
yields the most reasonable solution. The derived polar magnetic 
field strength is about twice the value previously published for 
the magnetic field mo dulus, a result that reappe ars for the MOST 
observations of y Equ dGruberbauer et al.l2008l) . The derived po- 
sition of 10 Aql deviates by less than 2<x from the luminosity 
based on the Hipparcos parallax. We conclude that this might be 
due to an inaccurate parallax measurement, but caution that our 
results must be considered with care as discussed in section [4] 
We also note that the spectroscopic data obtained simulaneously 
with MOST shows evidence for a pulsation node at a smaller op- 
tical depth (i.e., higher in the atmosph ere) than considered in the 
current models dSachkov et al.ll2008b . Therefore, the modeling 
effort and derived results presented in this work must be consid- 
ered as only a first attempt to understand the pulsation of this 
star. Still, our results presented here also reflect state of the art 
modelling efforts for pulsating magnetic stars and indicate the 
strong need for improved theoretical concepts. 

Despite the most complete and precise photometric cam- 
paign ever obtained for this star, the oscillation spectrum of 
10 Aql remains (compared to other, better studied members 
such as y Equ or HR 1217) poorly understood. Nevertheless, the 
MOST observations have certainly shed some new light on this 
star and provided new insights that might help to unravel the 
mysteries of this roAp star in future observations and modeling 
efforts. 

Acknowledgements. We are thankful to Mikhail Sachkov, Tanya Ryabchikova, 
Oleg Kochukhov and Nicole Nesvacil for discussion and collaboration with si- 
multaneous spectroscopy. It is also a pleasure to thank Chris Cameron for valu- 
able discussion on the theory of frequency spacings in roAp stars. DH, MG, 
MH, WWW, TK, PR and RK are supported by the Austrian Science Fund 
(FWF P17580). The Austrian participation in the MOST project is funded by 
the Austrian Research Promotion Agency (FFG). HS is supported by the 21st 
Century COE programme of MEXT, Japan. DBG, JMM, AFJM, and SR ac- 
knowledge funding from the Natural Sciences & Engineering Research Council 
(NSERC) Canada. This research has made use of the SIMBAD database, oper- 
ated at CDS, Strasbourg, France. 



Kochukhov, O., Landstreet, J. D., Ryabchikova, T., Weiss, W. W., & Kupka, F. 

2002, MNRAS, 337, LI 
Kurtz, D. W. 1978, Informational Bulletin on Variable Stars, 1436, 1 
Kurtz, D. W. 1982, MNRAS, 200, 807 

Lenz, P. & Breger, M. 2004, in IAU Symposium, Vol. 224, The A-Star Puzzle, 

ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 786-790 
Matthews, J. M., Kurtz, D. W., & Martinez, P. 1999, ApJ, 511, 422 
Preston, G. W. 1970, PASP, 82, 878 
Reegen, P. 2007, A&A, 467, 1353 

Reegen, P., Kallinger, T., Frast, D., et al. 2006, MNRAS, 367, 1417 

Rowe, J. F., Matthews, J. M., Kuschnig, R., et al. 2006, Memorie della Societa 

Astronomica Italiana, 77, 282 
Ryabchikova, T., Wade, G. A., Auriere, M., et al. 2005, A&A, 429, L55 
Ryabchikova, T. A., Savanov, I. S., Hatzes, A. P., Weiss, W. W., & Handler, G. 

2000, A&A, 357, 981 
Sachkov, M., Kochukhov, O., Ryabchikova, T., et al. 2008, Contributions of the 

Astronomical Observatory Skalnate Pleso, accepted; (arXiv:0712.1340) 
Saio, H. 2005, MNRAS, 360, 1022 
Saio, H. & Gautschy, A. 2004, MNRAS, 350, 485 
Tassoul, M. 1980, ApJS, 43, 469 

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023 



References 

Babcock, H. W. 1958, ApJS, 3, 141 

Baglin, A. & The COROT Team. 1998, in IAU Symposium, Vol. 185, New Eyes 
to See Inside the Sun and Stars, ed. F.-L. Deubner, J. Christensen-Dalsgaard, 
& D. Kurtz, 301-+ 

Balmforth, N. J., Cunha, M. S., Dolez, N., Gough, D. O., & Vauclair, S. 2001, 
MNRAS, 323, 362 

Belmonte, J. A., Kreidl, T. J., & Martinez Roger, C. 1992, Informational Bulletin 

on Variable Stars, 3752, 1 
Belmonte, J. A., Martinez Roger, C, & Roca Cortes, T. 1991, A&A, 248, 541 
Bigot, L. & Dziembowski, W. A. 2002, A&A, 391, 235 

Bigot, L. & Weiss, W. W. 2002, Communications in Asteroseismology, 141, 26 
Cunha, M. S. 2006, MNRAS, 365, 153 

Gruberbauer, M., Saio, H., Huber, D., et al. 2008, A&A, 480, 223 
Guenther, D. B. & Brown, K. I. T. 2004, ApJ, 600, 419 

Guenther, D. B., Demarque, P., Kim, Y.-C, & Pinsonneault, M. H. 1992, ApJ, 
387, 372 

Hatzes, A. P. & Mkrtichian, D. E. 2005, A&A, 430, 279 

Heller, C. H. & Kramer, K. S. 1990, MNRAS, 244, 372 

Huber, D. & Reegen, P. 2008, Communications in Asteroseismology, 152, 77 

Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 

Kallinger, T., Reegen, P., & Weiss, W. W. 2008, A&A, accepted; 

(arXiv:0801.0683vl) 
Kochukhov, O. 2004, ApJ, 615, L149 
Kochukhov, O. 2006, A&A, 446, 1051 
Kochukhov, O. & Bagnulo, S. 2006, A&A, 450, 763 



